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The most general amplitude for the radiative pion decay n — > ez/7 including terms beyond V- 
A theory is considered. The experimental constraints on the decay amplitude components are 
discussed. A model independent presentation of the results of high statistics and high resolution 
experiments is suggested. 
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Due to the strong helicity suppression of the n ev 
decay, the corresponding radiative decay 

TT — > eh'"/ (1) 

is regarded as a valuable source of pion structure infor- 
mation and a sensitive test anomalous interactions 0, |2| ■ 
Low's theorem Q allows the part of the tt ev^ am- 
plitude depending on the photon energy w as uj~^ and 
Lo^ to be unambiguously defined by the pion decay con- 
stant /tt = 131 MeV. This contribution to the decay 
amplitude, called inner bremsstrahlung [IB), is associ- 
ated with the tree diagrams shown in Figs. ^) and^^)- 
Generally, one can construct the remaining, structure de- 
pendent {SD) part of the amplitude (Fig. ^)) using the 
photon tensors J"^" = fc'^e'' - k^'t'' ± ie'"'"^fc"e'' where 
e'^ is the photon polarization vector and ± specify the 
photon helicities, the lepton currents L = Uv{\ + J^)ue, 
LA" = u^{l + j^)^t'ue, L'"' ^Uu{l+ -f^)a>"'ue, and three 
available 4-momenta k'^, q'' = + t'', and = — t^, 
where k^^, and are 4-momenta of photon, electron, 
and neutrino, respectively. Using the Dirac equation and 
properties of 7-matrices the full list of such amplitudes: 



TH^af^b^L, T^''a''b''LPcP, T'^ aPb"" LP" cP (F , J^^^a'^L", 
^ aP-hPL^P , and J^'^ L^'^ , may be reduced to four am- 



plitudes: 
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Tl''qP{qP LP" + 711,1") 



(2) 
(3) 

Tl"qPr"L (4) 



where me and m^r are the electron and pion masses, re- 
spectively. 




FIG. 1: The tt — > ewy diagrams 



With these ingredients, the full tt ~f ev^ amplitude 
may be presented as : 
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FtT + FsS\ (5) 



The first term in the sum is the IB amplitude. The 
structure dependent contributions are parameterized by 
form factors Fy , Fa, Ft, and Fs- In Eq. the vector 
and axial- vector form factors, Fy and Fa, are defined in 
accordance with the Particle Data Group (PDG)Q and 
the tensor form factor Ft uses the definition of Ref . |^ . 
In the most general case, form factors are functions of 
the kinematic variables x = 2fcg/m^ and A = kl/kq. A is 
expressed in terms of photon, x = lE^lniT^, and electron, 
y — 2Ee/mT^, center of mass energies asX={x + y— 1 — 
re)/x, where — TOg/m^. Neglecting the electron mass 
TOe, one can relate A to the e-7 opening angle in a simple 
form, A = j/sin2(6'ii;V2) = sin2(6'i;''V2), where e[^^ and 
9^^'^ are angles between e and 7 momenta in the pion and 
lepton pair centers of mass, respectively. In this paper 
we will assume that the form factors do not depend on 
A. Depending on the possible origin of the amplitudes, 
e.g. Tt"qPqPL''P cx (1 - x)T , T'^"qPkPL"P oc xT , the 
form factors may have a strong dependence on q^ or , 
equivalently, on x. 

To calculate the tt ev^ decay rate, the following 
sums over polarizations can be used: 

X^{l-X)\^ -TeX^X 

x^{l-x){\-Xf +reX^{l-X) 
x^X{l - X) 
x^{l~xfX{\-X) 
-rex2(l-A)[(l-x)(l 
V^2;2(l-A) 
^,x\l-X)[{l-x)X- 
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where px 
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plitude V±, 



PTV+ — PV+V- 

Ep„jA|V(8m4) and pxv 
In the limit rrig = 0, each am- 
T , S is characterized by definite photon 
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and electron helicities resulting the absence of the 
interference with each other. In this approximation, the 
decay rate, normalized to the branching ratio, may be 
presented as 



dxdX 



^ 1.43-10" 



(1 - .t)2 + 1 1 - A 



X 
\2 



+ |/y + /A|V(l-a;)A^ 
+ \fv~fA\^x^l-x){l-\y 

fT + fs{l~x)^]x{l-X) 
+ i f^ + fl{l-x)']x'X{l-X)\ (13) 



In Eq. (|13|l we have scaled the form factors 

fv,A,T,s = irnl/2meU)Fv,A,T,s ~ 14:6Fv^a,t,s (14) 

The first, IB, term is free of unknown parameters. The 
terms proportional to |/y±/yip are traditionally referred 
as SD^ terms. The term proportional to fr + /s(l — x)^ 
describes the interference of the amplitudes T and S with 
IB. 

In the Standard Model Ft = Fs = 0. The = {1 - 
x)'m'^ variation of the form factors Fv{q^) and FA{q^) 
is not expected to exceed a few percent over the total 
TT ev-f phase space fJjSH- The conservation of vector 
current (CVC) hypothesis relates the vector form factor 
Fy to the lifetime of the neutral pion 



\FviO)\ ^ 1 ^77) ^ 0.0259 ± 0.0005 (15) 

a \ 7rm_o 



[1| or equivalently fv ~ 3.78. As result, the tt — > evj 
decay amplitude contains only one unknown parameter 
7 = Fa/ Fy which has to be determined empirically. 

The TT —f ev^ experiments 9, 10, ,11, 12], which used 
the stopped pion technique, measured decay only in a 
limited phase space region where SD^ contribution is 
strongly dominant. By counting the tt — s- ev^ decays 
one could determine only |1 + 7^ with a small correction 
due to 1 1 — 7p. Some extension of the investigated phase 
space region allowed one to resolve the ambiguity 

in the sign of 7 and find a unique value: 



7 = 0.41 ±0.06 



(16) 



P|. This result is derived with the assumption that the 
value of the vector form factor in Eq. H15|) is valid. The 
choice of the sign of 7 was also confirmed by the 7r+ — > 
e^ve^e~ measurement jl3j . 

The ISTRA experiment [H Il5j| . in which pion de- 
cays were studied in flight, investigated the tt — )■ ev^ 
events over a much larger phase space region: 0.3 < 
X < 1, 0.2<A<1. The measured tt — > eu^ branch- 
ing ratio 1.61 ± 0.23, in units of 10"^, was found sig- 
nificantly smaller than the expected sum 2.41 ± 0.07 
of the IB (1.70), SD+ (0.67 ± 0.07), and SD- (0.04) 



contributions. The fit to the experimental distribution 
by the linear combination of the functions IB{x, A) and 
SD^{x,X) indicated a non-physical contribution of the 
SD-, BRsD- = -0.58 ± 0.20, to the total branching 
ratio. This result was interpreted Q as a possible indi- 
cation of a tensor amplitude with Ft = —0.0056 ±0.0017 
{/t = -0.82 ± 0.25). It was shown that this hy- 
pothesis not only does not contradict the results of the 
previous experiments fiol fill fl^ but even improves 
their consistency. 

Analysis of the new high statistics measurement of 
about 70,000 tt ey^ decays performed by the PI- 
BETA coflaboration hj is currently in progress. The 
preliminary results also indicate a deficit in the observed 
TT — > ev^ decays 18], however, a 2-parameter fit for 
Fa and Ft in the phase space region > 55.6 MeV 
and Ee > 20 MeV gives a much smaller value for 
Ft = -0.0017(1) (/t = -0.25 ±0.02) g?]. 

The PIBETA preliminary result calls for a more ac- 
curate study of the possible phenomenological explana- 
tion of the anomaly in the tt — > evj branching ratio. In 
this analysis we will assume that the discrepancy in the 
results are not caused by the experimental error nor by 
the statistical fluctuation. Since the measured number of 
events is smaller then expected, it can not be explained 
by some unknown background. Modification of the de- 
cay amplitude is required. To make an estimate, one can 
integrate Eq. H13|) over the ISTRA phase space region 
assuming constant form factors: 

10^BR^_,e,.7 = 1-70 + 0.327)1 + 7]^ + 0.169)1 - 7)^ 
+ 0.833/t + 0.212/1 + 0.099/s + 0.013/| (17) 

Depending on the value of 7, the total SD^ contribution 
to the decay rate can vary from 0.61 for 7 = 0.25 to 
0.79 for 7 = 0.52 assumptions about the form 

factors Fy and Fa, e.g. strong dependence on g^, vio- 
lation of CVC, big imaginary parts can not significantly 
change this estimate, because in Ref. they measured 
the structure dependent contribution to be 0.39 ± 0.07 
in the Ee > 56 MeV, E^ > 30 MeV, ei'^'^ > 132° region 
of the ISTRA total phase space. This value establishes 
a lower bound of 2.1 x 10"'' for the tt — > ei'j branching 
ratio if Ft = Fs = 0. Disagreement with the ISTRA 
result (1.61 ±0.23) x 10~^ clearly indicates the necessity 
to consider T and S contributions. 

From the quadratic dependence of the branching ratio 
on form factors one can derive the maximal reduction of 
the TT ev^ branching ratio which may be caused by /t 
or fs 



ABR: 
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= -0.19 
= -0.82 



{fs - -3.81) 
(/t = -1.95) 



(18) 
(19) 



One can see that the amplitude S with constant form fac- 
tor fs can not explain the ISTRA result. The possible de- 
pendencies of fs on X, e.g. fs = ff^x or fs = ff\l-x), 
does not change this conclusion. Therefore, the low value 
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TABLE I: The comparison of form factors frix) and ,fs{x) with the available experimental data. The form factors are 
parameterized by the constant /. For the combined analysis of the experiments ITfll ITU H^ . the estimate of / and the 
corresponding consistency of the data ^re shown. For ISTRA, indicate the agreement between the simulated and 
measured values of BRib and BRqq-. Predictions for the systematic shift of 7) A7 and fx, which one can obtain from the 
preliminary PIBETA fit are based on the estimation of the / in the simulation of the ISTRA fit. The predicted values of the 
/t may be compared with the preliminary experimental result /t « 0.25 0|. 



Considered form factors 
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-9.2l^,^ 
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0.0 


-0.01 
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q c+a.s 


3.35 


n 9+0.7 
^•^-0.8 


6.54 


0.01 


-0.29 ±0.07 



of the TT — > eiyj branching ratio, measured by ISTRA ex- 
periment implies the presence of the tensor amphtude T 
with a negative form factor. 

To evaluate a possible x-dependence of the form factor 
frix) we can simulate the available experimental data. 
For experiments H EIH El a combined has been 
constructed as a function of the "true" ratio of the axial 
to vector form factors 7 and a hypothetical form factor 
fT,s{x)- Minimization of the allows one to evaluate 
the significance of the /t,s for improving the consistency 
of the experimental data. The details of the analysis are 
described in Ref. _16]. 

To simulate the ISTRA results, the assumed IB and 
SD^ event distribution was modified by the hypothetical 
T (or S) contribution and then it was fit by a linear com- 
bination oi IB, SD^, and SD^ . The branching ratios of 
IB and SD~ were then compared with the ISTRA val- 
ues BRiB = 1.62 ± 0.20 and BRsd- = -0.58 ± 0.20 (in 
units of 10~^) using the function = {^js + '^sd- + 
2p6iBSsD-)/ii - P'), where Sib = [BRib - 1.62)/0.20, 
Sgjj- = {BRsD- + 0.58)/0.20 and p = 0.67 is the es- 
timate of the correlation factor. The dependence of 
on the studied form factor, allows us to simulate the re- 
sult of the determination of this form factor with ISTRA 
data. 

The form factors obtained in the ISTRA fit simulation 
were used to predict the results of the PIBETA fit in the 
kinematic region Ee > 20 MeV, > 55.6 MeV. 

In the simulation of the ISTRA and PIBETA fits, the 
detection efficiency was assumed to be constant in the 
phase space regions of interest. The stability of the re- 
sults against the possible efficiency dependence on x and 
A was specifically checked. 7 = 0.5 was assumed in both 
simulations. 

Several dependencies of /t and fs on x were tested. 
The correlated contribution of the form factors fx = 
/(I — x) and fs = f may be associated with amplitude 
(fc^e^ — k'^ e^)q^ L^'^ . The results are summarized in 



Table D 

One can see that form factor fs can not improve the 
consistency of the experiments U 0, 0, ^2 , does not 
fit well the ISTRA result, and does not effect the PI- 
BETA fit. In other words, there is no evidence of the 
structure dependent amplitude S in the available exper- 
imental data. 

The tensor form factors fx ~ f and fx — fx al- 
most perfectly simulate the ISTRA results. However they 
strongly disagree with the preliminary PIBETA fit. 

A form factor fx = /(I — x) is not in good consis- 
tency with the ISTRA data. Its value found in the IS- 
TRA simulation fairly well reproduces the total deficit 
of the TT — !■ evj decays. However, the distribution of 
the simulated missing events is better approximated by 
the IB{x,X) than by the SD~ distribution. Taking 
into account the correlation between the IB(x, A) and 
SD~{x, A) distributions, the low statistics of the ISTRA 
experiment (about 100 events), and the possible uncer- 
tainties in the simulation procedure, we should not rule 
out form factor fx — /(I — x) from the consideration. In 
addition, this form factor seems to be the only one which 
does not strongly contradicts both the ISTRA and the 
preliminary PIBETA results. 

We may also point out that the value of 7 which one 
can obtain in the PIBETA preliminary fit has only a weak 
dependence on the actual form factors fx{x) or fs{x). 

The value of the constant tensor form factor Fx — 
-0.0115 ± 0.0033 {fx = -1.7 ± 0.5) obtained in the dis- 
cussed simulation of the ISTRA fit is by factor of two 
larger in absolute value than the estimate in Ref. Q . In 
Ref. 0, = -0.0056 ± 0.0017 was estimated by com- 
paring the measured BRsd- with the calculated inter- 
ference between tensor and IB amplitudes. In this paper 
we took into account the terms quadratic in Fx and re- 
quired the part of total contribution of Fx which may 
be approximated by functions IB(x,X) and SD~ {x, X) 
to equal the measured sum of IB and SD~ contribution. 
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This promises to be a better approach. 

It is needless to eraphasize that the indirect analysis 
described above can not substitute the real fit to the ex- 
perimental data with the well determined detector effi- 
ciencies and resolution. One of the main goals of this pa- 
per is to stress that the analysis of the event distributions 
over A may provide model independent determination of 
the form factors. As one can see from Eq. each 
contribution to the decay rate have the same (known) A- 
distribution independent of the value of x and the actual 
value of the form factors. Since 1 — A = (1 — A)^ + A(l — A), 
only three values may be independently determined for 
each X 

a[x) = \fv + lA?x''{l^x) (20) 
h{x) = 2x [Jt + + Us + fh^){l - xf] (21) 
c{x) = X [\fv-fA?x''{l-x) + 4/t -I- 4/s(l-a;)2] (22) 

from the analysis of the events distribution over A, 

dN{x)/d\ cx a{x)\^ + 2h{x)\{l-X) + c{x){l-\f (23) 

We assume that the inner bremsstrahlung contribution 
IB{x){l - A)/(3A) with IB{x) = 3[(1 - xf + l]/x 13 
is properly subtracted from the experimental data. Ob- 
viously, the nonzero value of h{x) or negative value c{x) 
is evidence for the presence of the anomalous amplitudes 
T and/or S. Also, c(l) is directly related to the tensor 

form factor f^/"°^ = c(l)/4 at = 0. To illustrate 
the possible sensitivity of the method, the functions of 
a{x), b{x), c{x) are shown in Fig. |21 Functions b{x) and 
c{x) were calculated for three different tensor form factors 
obtained from the simulation of ISTRA fit: /t = — 1.7, 
/t = -2.8a;, and /r = -4.1(1 - x). Fig. |21shows these 
three cases are indistinguishable for x « 0.6. However, 
the measurement at a; > 0.7 allows one to make a choice. 

The measurement of a(x) allows one to determine the 
\fvix) + fA{x)\ as a function of x. Assuming that the 
vector form factor is known from the CVC hypothesis 
and chiral perturbation theory calculations ^ Q one can 
extract, in principal, fA{x), frix), and fs{x) from the 
measured values a{x), b{x), and c{x). 

In the above analysis of the tt — > ev^ decay rate, it 
was implicitly assumed that the form factors /t and fs 
are real and that there is no right handed neutrino in 
the decay amplitude. Such additional possibilities may 
be easily accounted for by the following modification of 
the form factors in Eqs. (|13I21I22|I : 



/t,5 
/t,5 



Here, are the form factors of the amplitudes with 
right-handed neutrino. Alternatively, /t,s and g may 
be considered as independent variables constrained only 
by the inequality f^g> /t,s x /t,s- For the tensor form 

factor at = 0, this inequality may be tested by 



ReifT,s) 
l^e(/T,s) 



\Im{f'_ 



T,S) 



(24) 
(25) 
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FIG. 2: The dependence of the functions a{x), b{x), c{x) on 
X, assuming fv — 3.78, /a ~ 0.5 fv, and fs — 0. For the 
functions b{x) and c{x), solid lines correspond to /t ~ —1-7, 
dashed lines to /t ~ —4.1(1 — x), and dashed-dotted line to 
/t = —2.8x. The inner bremsstrahlung contribution IB{x) is 
displayed for comparison. 



the measurement of 6(1) and c(l). With such an exten- 
sion of the decay amplitude we can no longer determine 
all form factors in a model independent way. However, 
the functions a{x)^ b(x), and c{x) still contain the com- 
plete experimental information which may be used in fur- 
ther theoretical analysis. In other words, the functions 
a{x)/x^{l — x), b{x)/2x, and c{x)/x provide the model 
independent representation of the results of the n — > ev"f 
experiment. Neglecting radiative corrections and correc- 
tions due to the charged lepton mass, me, this statement 
is valid in the general case and may be violated only by 
the hypothetical dependence of the form factors on on 
the e-7 opening angle. 

The discussion of the possible sources of the anomalous 
amplitudes S and T are beyond the scope this paper. 
Here, we will only mention that numerous attempts to 
explain the ISTRA result in the framework of the Stan- 
dard Model, e.g. due to the radiative corrections to the 
TT ev"f decay [22, were unsuccessful. 

In this paper the most general tt — > ev^ amplitude 
was considered. The model independent formulation of 
the results of the high statistics and high resolution ex- 
periments is suggested. If the final PIBETA result will 
validate the deficit of the tt — > ev^ decays observed by the 
ISTRA experiment, it will be a solid argument in favor of 
presence of the tensor component in the tt —t ev^ struc- 
ture dependent amplitude. Comparison of the available 
experimental data with the preliminary PIBETA results 
indicate a possible strong dependence of tensor form fac- 
tor on q^ . The large statistics accumulated by the PI- 
BETA experiment allows one to expect that the problem 
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of the TT — > ei/7 decay will be resolved soon. 
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